Orbital-assisted metal-insulator transition in VO2 
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We found direct experimental evidence for an orbital switching in the V 3d states across the 
metal-insulator transition in VO2. We have used soft-x-ray absorption spectroscopy at the V -t/2,3 
edges as a sensitive local probe, and have determined quantitatively the orbital polarizations. These 
results strongly suggest that, in going from the metallic to the insulating state, the orbital occupation 
changes in a manner that charge fluctuations and effective band widths are reduced, that the system 
becomes more 1-dimensional and more susceptible to a Peierls-like transition, and that the required 
massive orbital switching can only be made if the system is close to a Mott insulating regime. 

PACS numbers: 71.10.-w, 71.30.+h, 71.70.-d, 78.70.Dm 



The problem of metal-insulator transitions (MIT) in 
transition metal compounds attracts considerable atten- 
tion already for a long time. Among the best studied of 
such systems are the V oxides, especially V2O3 and VO2 
P, Q- The long-standing problem in these systems is 
the relative role of electron-lattice interactions and cor- 
responding structural distortions versus electron corre- 
lations. This problem is especially acute for the MIT 
in VO2, which was described either as predominantly a 
Peierls transition Q or as a Mott-Hubbard transition Q . 

An intriguing aspect that has largely been neglected 
in the discussions about MIT in TM oxides is the possi- 
ble role of magnetic correlations and especially the orbital 
structure of constituent TM ions || . Very recently, a the- 
oretical model for spinels such as MgTi 2 04 and CuIr 2 S4 
has been proposed in which a specific orbital occupa- 
tion effectively leads to the formation of one-dimensional 
bands, making the systems^in a natural manner, suscep- 
tible to a Peierls-like MIT |fj. 

The crystal and electronic structure of VO2 is in this 
respect more intricate. The MIT in VO2 is a structural 
transition from the high-temperature rutile (R) structure 
to a monoclinic (Mi) structure, in which there appears si- 
multaneous dimerization in each V chain along the c-axis 
and a twisting of V-V pairs due to an antiferroelectric 
shift of neighboring V atoms towards the apex oxygens 
(lying at the axis perpendicular to the crystal c-axis), see 
Fig. 1. As argued already long ago by Goodenough 0, 
one should discriminate between two types of orbitals and 
corresponding bands: dy-orbitals/bands, formed by the 
i2 S -orbitals with strong direct overlap with the neighbor- 
ing V in the chains, and 7r*-orbitals/bands, made of the 
two other t2g-orbitals. In the R-phase, the d\\ band over- 
laps with the 7r*-band, resulting in a orbitally isotropic 
metallic state, see Fig. 2. The twisting in the Mi -phase 




FIG. 1: Crystal structure of VO2 in the metallic rutile (R) 
phase and in the insulating monoclinic (Mi) phase. The ar- 
rows show the direction of the displacements of the V ions in 
the Mi-phase. The a, b and c-axes are defined with respect 
to the rutile structure. 



increases the V - apex O hybridization and moves the ir* 
band up, so that only the dn band is occupied. The later 
one then becomes split by the dimerization, leading to 
the insulating state. 

Many theoretical ab-initio studies were performed to 
test the Goodenough picture. LDA calculations indicated 
indeed that the dy band becomes more occupied in the 
Mi-phase 0, 0, but they failed to reproduce the insu- 
lating state. An L(S)DA+U approach predicted even 
more dramatic changes in the orbital occupations, but 
unfortunately it did not give the metallic solution for the 
R-phase |lfj . A change in the orbital occupation was 
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FIG. 2: Top panel: schematic electronic structure of VO2 
according to Goodenough Q; bottom panel: definitions of 
the relevant V 3d t^g orbitals used in this work, drawn in the 
(110) plane spanned by the a, b and c-axes of Fig. 1. Site 1 
and 2 are related by a 90° rotation around the c-axis. 



also obtained in the exact diagonalization study for a 
three-band Hubbard model using finite size clusters . 
Very recently, various LDA+DMFT methods have been 
applied to explain the MIT, and also here the orbital oc- 
cupations are an important issue [ToL IT3 . Il3l | . In view of 
the fact that the orbital occupation and changes thereof 
are central to the MIT theories for VO2, it is quite sur- 
prising that an experimental proof of it is lacking^ an 
aspect which is also relevant for other materials [lj, [l5j . 

In this paper we give a direct experimental evidence 
of this orbital redistribution at the MIT in VO2. We 
present a polarization-dependent x-ray absorption spec- 
troscopy (XAS) study on VO2 single crystals at the V 
£2,3 (2p — ► 3d) edges. Here we make use of the fact 
that the Coulomb interaction of the 2p core hole with the 
3d electrons is much larger than the 3d t^g band width, 
so that the absorption process is strongly excitonic and 
therefore well understood in terms of atomic-like tran- 
sitions to multiplet split final states subject to dipole 
selection rules. This makes the technique an extremely 
sensitive local probe 0, 0, 0] , ideal to study the or- 
bital character 0, 0] of the ground or initial state. For 
our experiment on VO2, we redefine the orbitals in terms 
of cr, 7r, or S with respect to the V chain as shown in Fig. 
2. The cr orbital is then equivalent to the du, and the w 
and 5 to the it*. The transition probability will strongly 
depend on which of the cr, tt, or 5 orbitals are occupied 
and on how the polarization vector E of the light is ori- 
ented. Our measurememts reveal a dramatic switching 
of the orbital occupation across the MIT, even more than 
in V2O3 0, indicating the crucial role of the orbitals and 
lattice in the correlated motion of the electrons. 

Single crystals of VO2 with Tmit = 67 °C have been 
grown by the vapor transport method |2£|. The XAS 
measurements were performed at the Dragon beamline 
of the NSRRC in Taiwan. The spectra were recorded us- 
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FIG. 3: Experimental V 1/2,3 XAS spectra of VO2 in the in- 
sulating Mi-phase (top panel, T=30°C) and metallic R-phase 
(bottom panel, T=100°C), taken with the light polarization 
E || c (solid lines) and E _L c (dashed lines). The metal- 
insulator transition temperature is 67 °C. 



ing the total electron yield method in a chamber with a 
base pressure of 3x10" 10 mbar. Clean sample areas were 
obtained by cleaving the crystals in-situ. The photon 
energy resolution at the V £2,3 edges (hv k, 510 — 530 
eV) was set at 0.15 eV, and the degree of linear polar- 
ization was sa 98%. The VO2 single crystal was mounted 
with the c-axis perpendicular to the Poynting vector of 
the light. By rotating the sample around this Poynting 
vector, the polarization of the electric field vector can be 
varied continuously from E \\ c to E _L c. This measure- 
ment geometry allows for an optical path of the incoming 
beam which is independent of the polarization, guaran- 
teeing a reliable comparison of the spectral line shapes 
as a function of polarization. 

Fig. 3 shows the V £2,3 XAS spectra of VO2 taken in 
the insulating Mi-phase (top panel), and in the metallic 
R-phase (bottom panel). The general lineshape of the 
spectra is similar to that in earlier works, each of which 
only reports spectra for one particular polarization |2ll 
I22] ]. In our work, we have measured for each phase the 
spectra with two different light polarizations, namely E \\ 
c (solid lines) and file (dashed lines). We observe a 
clear polarization dependence for the insulating phase. 
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FIG. 4: Theoretical simulations for the polarization dependent V L2,3 XAS spectra. Left panel: best fit to the experimental 
spectra, using orbital occupations as indicated in Table I. Right panel: simulations for the V 3d 1 ion assuming a pure a, tv, or 
S orbital occupation. Top part of the panels: simulations for the insulating Mi-phase. Bottom part of the panels: idem for the 
metallic R-phase. 



By contrast, the polarization dependence is quite weak 
for the metallic phase. Fig. 3 shows for each phase also 
the dichroic spectrum, i.e. the difference between the 
spectra taken with the two polarizations. One now can 
see that the dichroic spectrum of the insulating phase 
has not only a larger amplitude, but also a very different 
lincshape than that of the metallic phase. 

To extract information concerning the orbital occupa- 
tion from the V £2,3 XAS spectra, we have performed 
simulations of the atomic-like 2p 6 3d 1 — > 2p 5 3d 2 transi- 
tions using the well-proven configuration interaction clus- 
ter model [IE [HO. The method uses a V0 6 cluster 
which includes the full atomic multiplet theory and the 
local effects of the solid. It accounts for the intra-atomic 
3d-3d and 2p-3d Coulomb interactions, the atomic 2p and 
3d spin-orbit couplings, the O 2p - V 3d hybridization, 
and local crystal field parameters. Parameters for the 
multipole part of the Coulomb interactions were given by 
the Hartree-Fock values [lif. while the monopole parts 
(Udd, U p d) as well as the O 2p - V 3d charge transfer- 
energy were determined from photoemission experiments 
|23| . The one-electron parameters such as the O 2p - V 3d 
and O 2p - O 2p transfer integrals as well as the local crys- 
tal fields were extracted from the LDA band structure 
results for the crystal structure corresponding to each 
of the two phases of VO2. The simulations have been 
carried out using the XTLS 8.0 programm [IH Hil Eij . 

Fig. 4 shows the results of our theoretical simulations 
of the spectra. In the top part of the right panel we have 
simulated the insulating Mi-phase spectra for the follow- 
ing three scenarios: the V 3d 1 ion is set either in the pure 
<7, 7r, or 5-orbital symmetry. One can clearly observe that 
the different orbital symmetries will lead to very differ- 
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TABLE I: Symmetry and orbital occupation of the 3d shell 
of VO2 in the Mi- and R-phase. 



ent spectra with quite different polarization dependence. 
One can notice that the cr-orbital scenario resembles the 
experimental spectra the most, especially when one fo- 
cuses on the most excitonic part of the spectrum, namely 
between 512 and 516 eV. In a simulation with the V ion 
in the ground state that belongs to the true local sym- 
metry of the Mi-phase, we find even a better fit to the 
experimental data as shown in the left panel. The corre- 
sponding orbital symmetry, as listed in the left column 
of Table I, has indeed overwhelmingly the a character 
(0.81), and only very little tt (0.10) and S (0.09). 

We have also simulated the spectra for the metallic 
R-phase, again for the three scenarios in which the V 
ion set to have either the pure er, 71", or (5-orbital sym- 
metry. The bottom part of the right panel shows that 
each scenario results in quite different spectra and po- 
larization dependence. We also note that each of the 
R-phase scenario gives spectra different from the corre- 
sponding Mi -phase, simply because of the differences in 
the local electronic structure, resulting from the differ- 
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ent crystal structure. Important is now that none of the 
three scenarios of the R-phase give good agreement with 
the experimental spectra. Apparently, the V ion has an 
orbital symmetry which is very far from a pure a, or it, or 
S. We now approximate the initial state symmetry of the 
V ion by a linear combination of those three symmetries, 
and optimize the relative weights to obtain the best fit to 
the experiment, with the emphasis on the excitonic part. 
Fig. 4 (bottom left panel) shows that this state is built 
up of 0.33 a, 0.16 ir and 0.51 6 symmetries, see also the 
5th column of Table I. It seems thus that in the metallic 
phase the V orbital occupation is almost isotropic. 

In Table I we have also listed the 3d orbital occupation 
as found from the simulations of the experimental spec- 
tra. These numbers are not identical to the symmetry 
occupation numbers because of the covalency, i.e. the 
hybridization of the V 3d with the surrounding O 2p lig- 
ands. We now can compare our findings directly with the 
numbers from our LDA and LSDA+U calculations ||. 
We note that our LDA band structure is quite similar 
to the one published earlier 0, || , and that the occupa- 
tion numbers of our LSDA+U is in close agreement with 
the one published very recently [inj . For the insulating 
Mi-phase, we find that the orbital occupation, which is 
highly (cr) polarized, is well reproduced by the LSDA+U 
model but not so by the standard LDA, see Table I. On 
the other hand, for the metallic R-phase, we observe that 
the almost isotropic orbital occupation as experimentally 
determined is well reproduced by the LDA, but not so by 
the LSDA+U. It seems that the LDA tends to underes- 
timate the orbital polarization, which makes the method 
less suitable for the insulating phase. The LSDA+U, on 
the other hand, tend to overestimate it, which puts this 
approach in disadvantage for the metallic phase. These 
problems are likely to be related to the fact that the LDA 
cannot reproduce the insulating state in the Mi-phase, 
while the LSDA+U does not give the metallic state for 
the R-phase. Nevertheless, the general trend that the or- 
bital occupation is more cr-polarized in the Mi-phase is 
predicted correctly in both approaches. 

In comparing our experimental results with the DMFT 
calculations, we note that in one implementation of the 
standard LDA+DMFT method the change in orbital po- 
larization is too small which has been attributed to the 
fact that the insulating phase cannot be reproduced llCl . 
Very exciting is that an LDA+cluster /DMFT study 
has been very successful in reproducing the strong change 
in orbital polarization, indicating the importance of the 
/c-dependence of the self-energy correction. 

The significant outcome of our experiments is that the 
orbital occupation changes from almost isotropic in the 
metallic phase to the almost completely cr-polarized in 
the insulating phase, in close agreement with the three 
band Hubbard model 0] . This very strong change leads 
to a dramatic modification of the intersite exchange in- 
teractions with large consequences for the effective Hub- 



bard U for nearest neighbor charge fluctuations 0, Hl| 
and effective band widths Moreover, the orbital po- 
larization change is such that VO2 in terms of its elec- 
tronic structure is transformed from a 3-dimcnsional to 
effectively a 1-dimensional system |(|. The V ions in 
the chain along the c-axis are then very susceptible to a 
Peierls transition. In this respect, the MIT in VO2 indeed 
has many features of a Peierls transition However, 
to get the dramatic change of the orbital occupation one 
also need the condition that strong electron correlations 
bring this narrow band system close to the Mott regime 
[4j. The MIT in VO2 may therefore be labelled as an 
orbital assisted "collaborative" Mott-Peierls transition. 

To conclude, we have found direct experimental evi- 
dence for an orbital switching in the V 3d states across 
the metal-insulator transition in VO2. The orbital oc- 
cupation in the insulating state is such that the effec- 
tive band widths are reduced and the system electroni- 
cally more 1-dimensional and thus more susceptible to a 
Peierls-like transition. 
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